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As a candidate scheme for controllably coupled qubits, we consider two quantum dots, each doped 
with a single electron. The spin of the electron defines our qubit basis and trion states can be created 
by using polarized light; we show that the form of the excited trion depends on the state of the qubit. 
By using the Luttinger-Kohn Hamiltonian we calculate the form of these trion states in the presence 
of light-heavy hole mixing, and show that they can interact through both the Forster transfer and 
static dipole-dipole interactions. Finally, we demonstrate that by using chirped laser pulses, it is 
possible to perform a two-qubit gate in this system by adiabatically following the eigenstates as a 
function of laser detuning. These gates are robust in that they operate with any realistic degree of 
hole mixing, and for either type of trion-trion coupling. 
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I. INTRODUCTION 



In the quest for a solid state quantum information pro- 
cessor, there is great attraction in combining the rela- 
tively long coherence times of spins with the speed and 
versatility of optical manipulation. There have been 
many papers describing different ways to embody a qubit 
by using the two levels of a confined spin-1/2 elec- 
tori > i i 2 4 3 j 4 i 5 4 6 The direct interaction between two such 
spin qubits is often quite weak, but they can be enhanced 
by exploiting degrees of freedom which lie outside the 
computation Hilbert space. For example, spin informa- 
tion can be transferred to spatial degrees of freedom in 
a double quantum dot structure^ or to photons in an 
electromagnetic cavityi£ Recent proposals&2*i2i have de- 
scribed ways in which spins might be coupled by using 
polarized light to selectively create trion states. Such 
spin selective charge excitations benefit directly from 
the recent progress in ultrafast optoelectronics: both 
the coherent manipulation of excitons in quantum dots 
(QDs) 1112 ' 1 - and spin selective optical transitions 1 ^ have 
been demonstrated. Moreover, once a gate operation 
is complete, it is possible to arrange that all popula- 
tion returns to the qubit subspace, and the quantum 
device therefore benefits from the robust coherence prop- 
erties 15,16,17 of the electron spin. 

In Ref. IToL we demonstrated that an entangling 
CPHASE gate can be performed between single spins on 
each of two adjacent QDs by the spin selective excitation 
of a single, delocalized, exciton. The required delocal- 
ization occurs when the dots are (near) resonant, and if 
they interact through the Forster energy transfer mech- 



anism 



18.19 



The proposal is based on the Pauli blocking 
mechanism 2 ^ and is valid only in the special case of no 
light-heavy hole mixing, which, though sometimes appli- 
cable^ is not generally true for most QD systems^ We 
shall here demonstrate an alternative method of perform- 



ing a two-qubit gate, which solves this problem. Follow- 
ing Ref. !). it is based on using chirped laser pulses to 
perform adiabatic quantum gate operations. 

To develop this model, the following steps are nec- 
essary. We shall first describe how the QDs employed 
in our quantum information implementation scheme are 
modeled. We describe the hole sub-band mixing in terms 
of a four band Luttinger-Kohn model^ and show how 
this affects the coupling of confined charge carriers to 
a laser field (Section [nj . We then consider a coupled 
QD system and present the dependence of the Forster 
transfer operator on the angular momenta of the exci- 
tonic states which it connects. We shall derive the form 
of the Forster transfer interaction for two coupled tri- 
ons and hence write down an effective Hamiltonian for 
two interacting QDs coupled to a laser field, in the pres- 
ence of hole mixing (Section II I II) . The proposal for per- 
forming quantum gates by creating trions adiabatically is 
then discussed (Section II VII . and we describe two differ- 
ent modes in which the two qubit gate can be operated. 
We shall demonstrate that the adiabatic scheme circum- 
vents the hole mixing problem and briefly discuss why it 
should also reduce phonon decoherence. State measure- 
ment, preparation and scalability will be discussed next 
(Section 01 and then we summarize (Section lv"T|l . 



II. SINGLE QUANTUM DOT MODEL 

Let us consider self-assembled QDs with strong con- 
finement along the growth direction z, which is also 
the QD symmetry axis. This type of QD can be pro- 
duced in materials such as InGaAs by using the Stranski- 
Krastanow method, 24 which may allow the realization of 
a controllably coupled many dot system. Such QDs ex- 
ist in the strong confinement regime, in which the typi- 
cal size, L, of the QD in the growth direction is of the 
order of 10 — 20 nm. In this regime, the Coulomb in- 
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teraction between charge carriers scales as 1/L, but the 
single-particle excitation energy has a 1/L 2 dependence. 
Excitonic wavefunctions can therefore be modelled by 
products of single particle electron and hole states, with 
Coulomb effects being introduced by using first order per- 
turbation theory. This approximation results in a shift 
to the excitonic energy but does not lead to the entan- 
glement of electron and hole. The effective mass and 
envelope function approximations reveal that the elec- 
tronic states inside the QD exhibit atomic-like symme- 
tries, which have been identified experimentally 2 ^ 

The wavefunction for a single particle in a QD can 
be described by a product of a Block function U, which 
has the periodicity of the atomic lattice, and an envelope 
function (f>, which describes the amplitude modulation 
of the wavefunction that is imposed by the confinement 
potential. Henceforth we shall only consider the lowest 
energy envelope function for both the conduction and 
valence bands (which has no nodes in both cases), and 
neglect any mixing with higher envelopes. This approx- 
imation is discussed in Ref. l26t we use it for clarity and 
our scheme does not depend on it; a more thorough de- 
sciption of the electronic structure of self-assembled QDs 
is presented in Ref. |U 

The eigenstates of the angular momentum operators, 
J and J z for the six hole states closest to the top of the 
valence band can be represented by^ 



The electron states are simply: 
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We have labeled the Bloch functions U by using the no- 
tation | J/j, J z ). The first two states correspond to heavy 
holes (hh); the next two are light holes (Ih) and the 
last two are split-off holes (so). The functions fi de- 
scribe the radial dependence of each Bloch function type 
i £ {hh, Ih, so}; a and /3 are the up and down spin states 
respectively. The X, Y, and Z represent orbital wave- 
functions as follows: 



(r\X) 
(r\Y) 
(r\Z) 



cos 9. 



(7) 
(8) 
(9) 



|l/2 e , 1/2) = g(r)\Sa), 
|l/2 e ,-l/2) = g(r)\S(3). 



(10) 
(11) 



(r| S) — 1/V47T is the isotropic orbital function and g(r) 
is the radial dependence of the electron's wavefunction. 
We now note that: 

[r\X)x(r\S)dil = J (r| Y)y (r| S)dQ = 



(t\Z)z(t\ S)dfl 



V3 



(12) 



where dJ7 is the infinitesimal solid angle. We also see 
that: 

J {r\X)y(r\S)dn = J (r| X)z (r| S)dn = (13) 
(r|Y)a;(r|5)tff2 = J (r| Y)z (r| S)dtt = (14) 
(r| Z)x (r| S)dil = [ (r| Z)y (r| S)dil = 0. (15) 



These relations will be important in the following discus- 
sion. 



A. Hole Mixing 

Most semiconductors exhibit mixing of the heavy 
and light hole sub-bands, which we shall now describe 
by using the Luttinger-Kohn model. 23 The electron 
eigenstates of bulk semiconductors may be characterized 
by the crystal momentum wave- vector, k = {k x , k y ,k z }. 
The coupling between light and heavy holes is described 
by a four band Luttinger-Kohn Hamiltonian^ so long 
as the split off holes are energetically distant enough that 
coupling to this band can be neglected^ In the basis 
{\J Z - +3/2) , | J* = +1/2) , \J Z = -1/2) , \J Z = -3/2)}, 
the Hamiltonian is written: 
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The uncoupled heavy hole Hamiltonian, Hhh is: 

h 2 k 2 z , s h 2 {k 2 x + kl) i 

^ = ^( 7l - 2 72 ) + ^-^( 7l + 72 ). 
The uncoupled light hole Hamiltonian, Hih is: 

The mixing parameters are: 
V3fr 2 



(16) 
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\l2(k x - k ) - 2ij 3 k x k y ] 
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and 



iky). 



(20) 



71, 72 and 73 are the Luttinger parameters***** and mo 
is the free electron mass. 

For the quantum confined states that are found in 
semiconductor nanostructures, the crystal momentum is 
no longer a good quantum number. It is therefore neces- 
sary to replace each component of k with its expectation 
value, taken over the hole envelope function, <fi(r)£& That 



(r)V0(r)rfr = (k) 



(21) 



If the envelope function has a well-defined parity then 
(k x ) = (k y ) = (k z ) = 0, and so, from Eq. EDI b = 0. 
Therefore, the Hamiltonian, Eq. ^| decouples and acts 
in two separate two-dimensional Hilbert spaces. We find 
that the hole eigenstates are given by the two pairs: 



\h+) = y/l - e 2 \J Z = +3/2) +e\J z = -1/2) , 
\h'_) = y/l - e 2 \J Z = -1/2) -e\J z = +3/2) , (22) 



and 



\h-) 



\J Z 



-3/2) + e\J z = +1/2) , 



\ti + ) = y/l - e 2 \J Z = +1/2) -e\J z = -3/2) , (23) 

where e characterizes the degree of mixing. \h+) and 
\h—) are degenerate states that are predominantly heavy- 
hole like, and which are split from the second degenerate 
pair, and \ h'_) 1 which have predominantly light-hole 
character. The nature of mixing here is quite different 
to the type considered in Ref. |9j, where mixing between 
the \ J Z = +3/2) and \J Z = +1/2) (or \J Z = -3/2) and 
\J Z — —1/2)) states was assumed. 

To estimate the hole mixing parameter let us now 
consider a specific and very simple model: that of a 
hole bound in a parabolic potential in all three dimen- 
sions. 30 This external potential is defined by V(x, y, z) = 
(u> x x 2 + iJ^y 1 + w 2 z 2 )/7 1? 2& where ujj is the frequency 
of the trapping potential in the j = {x, y, z} direction. 
We shall make the axial approximation - i.e. we shall 
ignore terms which are not axially symmetric about the 
z-axis; this is a good approximation for GaAs^i (cor- 
rections to this approximation are discussed in Ref. |32J). 
Then the Luttinger-Kohn Hamiltonian fEq. H6|) becomes, 
in the {\J Z = +3/2) , \ J Z = —1/2)} basis (or equivalently 
the {\J Z = -3/2) , \J Z = +1/2)} basis): 
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W 

2uj t - AS,, 



(24) 



where, in order to simplify notation, we have defined: 



\/3( 72+73) 



uo x + ujy + uj z , AE h = (ujt - 3Wz 



and W = 

71 ' 



(UJ X -U)y), 



Assuming that the difference between the diagonal el- 
ements is much greater than the magnitude of the off- 
diagonal elements - i.e. 2AEh ^> W, then t« 1 and is 
given by: 



W _ \/3~7i(72 + 73) (Ux ~ u y ) 
2AE h ~ 472 {uj t - 3w z ) 



(25) 



As expected, the mixing is proportional to sub-band cou- 
pling and inversely proportional to the energy difference 
between the heavy and light hole sub-bands. Taking 
typical values for the Luttinger parameters for GaAs 
(71 = 6.8, 72 = 2.1 and 73 = 2.9), and estimating the 
different trapping frequencies for an anisotropic QD to be 
hjj x =10 meV, ftw y =11 meV and Hlo z =45 meV, leads 
to e w 0.1. We shall use this value throughout the rest 
of the paper. 



B. Relevant Trion States 

We shall now assume that the hole states \h'_/) and 
\h'_) are energetically distant enough from \h + ) and \h-) 
that they may be ignored in our calculations of quantum 
dynamics (in the parabolic well model considered in the 
previous section this splitting is of order 1 1 meV) . We are 
interested in using single excess spins in QDs to embody 
our qubit, and in exploiting spin-dependent exciton cre- 
ation to couple together spins in adjacent dots. A single 
spin becomes a trion stated 3 , following the creation of an 
exciton, and by using Eas. and 1231 we find that the 
trion eigenstates are: 



l^+> = !<%><B> \h+), 
\x-) = \S n ) ® \h-) , 



(26) 
(27) 



where Sn denotes two electrons in opposite spin states. 

We emphasize that the states \h+) and \h—) are not 
eigenstates of the J z operator - and this has profound 
consequences for both the coupling of charge carriers to 
the laser field as well as the Forster transfer interaction. 
We now specifically calculate the form of these two inter- 
actions. 



C. Interaction with a laser field 

The QD-light interaction for a (classical) laser pulse 
of amplitude E(t) and central frequency wj,(t) impinging 
on a single QD may be expressed in the dipole approx- 
imation 2 ^ through the following Hamiltonian operator: 



Ti L (t) = eE(t)r ■ hcos[ui L (t)t] 



(28) 



where f is the dipole operator, n is the polarization vec- 
tor of the light field, and e is the electronic charge. The 
time dependence of E and lj allows us later to introduce 
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chirped laser pulse shapes: we assume that the time de- 
pendence is slow compared with the oscillation period of 
the laser. 

We assume that the laser pulse has a spectral width 
which is narrower the typical QD level spacing. If we 
choose the laser frequency to be close to resonance with 
the ground state exciton, we can then consider an ide- 
alized model in which the dynamics is restricted only to 
the two qubit states defined as 



|0) = |-l/2 e 
11} = |l/2e> ; 



(29) 



and to the trionic states |x+) and \x—). 

By choosing the laser pulse to be directed along z and 
<j circularly polarized, we find that f • n = (x + iy). 
Eqs. HH HJ3 El and EHl then allow us to calculate the 
form of the interaction between the qubit electron spin 
and trion states. Let us first define the length 



(30) 



li = / fi(r)r 6 g(r)dr 



with i £ {lh, hh}. We also use the modified mixing pa- 
rameter e = e l ih anc j tfiQ j^a^j frequency 

IhhvS 



Q(t) 



2eE(t)l hh 
V6 ' 



(31) 



Then we find, for a quantum dot labelled by k: 



t (t) = n(t) cos[uj L (t)t} (|1) K (x+\ + e\0) K 



H.c.) 

In the absence of mixing (e = 0) only the spin up qubit 
state |1) is coupled to the laser field; the spin down state 
|0) is completely decoupled. However, for a finite amount 
of hole mixing the spin-selectivity of trion excitations (re- 
quired for the scheme of Ref . llOl) is no longer maintained 
(see Fig.nj and we must therefore consider an alternative 
gating strategy. 



III. COUPLED QUANTUM DOT STRUCTURE 

There are two principal interactions between trion 
states in adjacent quantum dots. The first of these is the 
direct Coulomb binding energy Vxx between two trions, 
which leads to the biexcitonic shift >2i The other interac- 
tion is the off-diagonal Forster coupling, which induces 
the transfer of electron hole pairs, through virtual pho- 
tons, and we discuss this in the following sections. 



+ 



vAAAA-* 



Mixing 




FIG. 1: Hole states in the presence of hole sub-band mix- 
ing and its effect on Pauli blocking. The arrows denote the 
allowed optical inter-band transitions for an incoming a + po- 
larized laser pulse. The left part of the figure represents the 
situation where no hole mixing is present and the right hand 
side shows the effects of mixing. 



found that the magnitude of the interaction is given by 



4-Treoe 



jpWtWi «ri).(r 2 ) 



- — « ri ).R)«r 2 >.R)){33) 



where R is the vector connecting the centers of the two 
QDs. The term (r^) represents the interband expectation 
value of the atomic position operator for dot i: 



U e (v)rU h (r)dr, 



(34) 



cell, dot i 



where the U e and Ut represent the Bloch functions for 
electrons and holes within dot i. The Wi represent the 
overlap of the envelope functions <p e and (j>h for dot i: 



IF 



4(r)^(r)dr. 



(35) 



An experiment would typically be performed on a pair 
of vertically stacked quantum dots, where the vector R 
lies along the growth direction z: by using this in Eg. 1331 
it is easy to show that the Forster interaction conserves 
the angular momentum of transferring excitons. How- 
ever, the nature of the exciton state before and after the 
transfer affects the magnitude Vf. By substituting the 
Bloch functions of Eqs. to into Eg. 1331 and using the 
relations of Eqs. ^| to El we can obtain the strength of 
the Forster coupling for excitons which are composed of 
electrons and holes of varying angular momentum. We 
defined 



Angular momentum dependence of the Forster 
interaction 



The general form of the Forster coupling Hamiltonian 
in QDs was discussed in detail in Ref. l35l where it was 



Mi. 



12ne e r R 3 



WiWuUlj, 



(36) 



for i,j £ {lh,hh}. Table I then shows the matrix ele- 
ment for all transitions which are induced by the Forster 
interaction. 
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State 1 {J z ,h, Jz, e ) 


State 2 (J z>h , Jz,e) 


Matrix Element 


Net J z 


3 . 1 

2 ' 2 
3 . 1 

2 ' 2 

1 . 1 

I. 5 ' ? 

2 'l. ? 

2 ' 2 

r. i 

2 ' 2 

r. i 

2 ' 2 

3 . 1 
2 ' 2 

3 . 1 

2' 2 


3. 1 

2 ' 2 
3 . 1 
2 ' 2 

1 . 1 

i 2 ' ! 
?'. l 

2 ' 2 

1 . 1 

"?:"i 2 

2 ' 2 

1 . 1 

2 ' 2 
1 . 1 

2' 2 


Mhh.hh 
Mhh,hh 

-4M !MA /3 
-4M !Mfc /3 

M lhl i h /3 
M lh}ih /3 
Mi h , hh /V3 
Mi h)hh /y/% 


-1 <-» -1 
+1 ^ +1 
0^0 
0^0 
0^0 

-1 <-> -1 
+1 ^ +1 
-1 <-> -1 
+1 <-> +1 



TABLE I: Relative size of the matrix element for exciton states coupled by the Forster interaction. M is defined in Eg. 1361 



B. Forster interaction for coupled trions 



Full Hamiltonian 



We next derive the form of the Forster transfer opera- 
tor T for two coupled trions. The most general definition 
of f is: 



J 



H.c. 



f= J2 v F J,k,: 

(37) 

where Vp ' is the size of the Forster matrix element 
which connects an electron and hole on dot 2 (whose an- 
gular momentum states are labelled with indices k and 
I) to an electron and hole and dot 1 (whose angular mo- 
mentum states are labelled with indices i and j). The 
only non-zero values of Vfj k l are given in Table 1, and 

we use that table to work out the effect of T on states in 
our two-dot system. For example, we have that: 



T\lx A 



= T 



— Mhh.hh 

eMihAh 

3 

±eMi K i h 
3 

We can now see that 



1 



3 Q 



1 9 - + 1 

"2V % ' + 2, 

+ 1 <f • 1 



(x+l|T|la:+) = M hh<hh + 

I n\rh\i \ ^ 2 Mlh,lh 

(x-Q\l \lx+) — - ■ 



Hh.lh 



(38) 

(39) 
(40) 



Similar calculations allow us to find all of the Forster 
coupling terms, which may be expressed by the following 
Hamiltonian, correct to first order in e: 

H F = M hhM {\Qx-) (z_0| + \lx + ) 
2M hh , h e 



V3 



(|lx_) (x+0| + (0x+\) + H.c 

(41) 



The total Hamiltonian of the two QD system in the 
presence of laser excitation may now be written as: 

K—a 7 b 

+ Yl V xx {\ V n){v^\). (42) 

The state |00) sets the zero of energy - and then the 
first term describes the Zeeman energy splitting 5 of the 
spin qubits; the second term represents the trion creation 
energy w Xtt where P Xk = \x+) K {x + \ + \x-) K (ar_| is the 
projection operator onto the single trion state located 
in the QD labeled by k; the next two terms are defined 
by Eqs. and ^2 respectively, and the last term de- 
scribes the static dipole-dipole binding energy between 
trions V XX 2& 

We now move to a frame which is rotating at the laser 
frequency ujl for both spin-trion transitions, and make 
the Rotating Wave Approximation (RWA). We then find 
that: 



H T (t) = ]T [*|l> B (l|+A«(t)P X( .+W^,„(i) 



Y v xx Wii) 



-TLf 
(43) 



where u> Xk has now been replaced by the time dependent 
detuning of the laser from the spin to trion transition 
energy A K (t) = uj Xk — cjr,(i). The charge-laser field cou- 
pling is now given by 



(\l) K (x + \+e\0) K (x-\+H.c). (44) 



The Hamiltonian, Eq. Ej3 spans a 16-dimensional 
Hilbert space. It is composed of the four computational 
basis states (|00), |01), 1 10) and 1 11}), eight single trion 
states, and four double trion states. In order to get a lit- 
tle more insight into the behavior of the system, we plot 
its eigenenergies as a function of the ratio A/fl in Fig. El 
There are three main groups of curves that are well sep- 
arated away from A/f2 = (i.e. when |A/f2| 3> 1); these 
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Single tripns 




-6 -4 -2 



FIG. 2: Eigenstate energy spectrum. The parameters are set 
as follows: 5/Q = 1, M hh , hh /Q, = Mi h , hh /Q, = 0.5, Vxx/£l = 
2, e = 0.1. 



groups are simply the four computational basis states, 
eight single trion states and four double trion states. As 
A/f2 approaches zero (i.e. when the laser becomes reso- 
nant with the spin to trion transition energies), the eigen- 
states become superpositions involving different numbers 
of trions, and many anticrossings can be seen in the eigen- 
state spectrum. 



IV. ADIABATIC QUANTUM GATES 

In this section, we shall describe how an adiabatic 
change in the ratio A/fi can allow us to follow the eigen- 
state curves. Not all of the computational basis states 
mix in the same way with trions when A/fl is small. 
Following Ref. |3, it is possible to use chirped laser pulses 
to slowly vary A/fi, which causes a non-trivial two qubit 
operation via trion state anticrossings. 



A. Chirped pulses 

Let us consider such a chirped laser pulse, where the 
detuning varies in time as follows: 



A(t) = -A (l-le-W^y 



(45) 



Ao represents the maximum detuning and ta is a param- 
eter characterizing the time variation of the detuning. 

The time dependence of the laser intensity is assumed 
to cause the Rabi frequency to take the following form: 



fl(t) = fl e 



(46) 



where, similarly, Qq represents the maximum Rabi fre- 
quency and Tn is a parameter characterizing its time 
variation. For a review of experimental methods for pulse 
shaping see Ref. 

We calculate the quantum dynamics caused by our 
Hamiltonian with these time-varying laser pulses by us- 
ing a numerical Schrodinger equation solver. If the adia- 
batic approximation holds true, population will return to 
the initial state at the end of the operation (non-adiabatic 
corrections limit the gate fidelity and are described be- 
low). However, the phase accumulated during the gate 
varies depending on the initial state. We characterize the 
gate by looking at the relative phase 9 gained when the 
pulse is applied to each of the four computational basis 
states in turn: 9 



no 



nx 



(47) 



where <f> n is the phase change of \n) during the gate oper- 
ation. The relative phase 9 is the part of the phase which 
is invariant under single qubit operations (see Ref. 37 for 
a detailed discussion) . We can see this by considering the 
effect of the following gates. First, 



Ui 



"'000 







-i<t>io 



(48) 



is performed on qubit 1 (in the |0), |1) basis), and then 

^2 = ^ q gi^OD-^lo) ^ ( 49 ) 

is performed on qubit 2. These two single qubit gate 
operations remove any phase picked up on the states |00), 
1 01} and 1 10), with |11) undergoing a net phase change 
of 9. In the basis |00), |01), |10), |11), a CPHASE gate 



Uc PHASE = 



10 

10 
10 
-1 



(50) 



This can be constructed from any gate operation in which 
9 = 7r, together with appropriate single qubit gates^i 

A strength of our approach is that a phase gate can 
be performed whether the interdot coupling takes the di- 
agonal form (biexcitonic interaction) or an off-diagonal 
form (Forster transfer). We can illustrate this by per- 
forming numerical simulations in two limits. First, we 
consider the case where Vxx > Mhh,hh = Mih,hh- In this 
case the required phase is given predominantly through 
the biexcitonic coupling. If we take the same parame- 
ters as used for Fig. [21 with chirped pulses characterized 



7 



by Tn = 3.55 ps, ta = 2.55 ps, A = 4.5 meV and 
fio = 8 mcV, then we can indeed obtain 9 — n. This is 
shown in Fig.0 which displays the time dependence of 9. 
The population of each computational basis state (follow- 
ing initialization into that state) is displayed in Fig. 
the chirped pulse moves some of the population out of 
each state (and into trion states) and then back again. 
The effect of the pulse is the same for 1 01) and 1 10) , as 
would be expected on further inspection of our Hamil- 
tonian. However, there is a distinct difference between 
the behavior of these two states and that of the states 
1 00) and 1 1 1) - and it is this difference which allows the 
relative entangling phase 9 to be picked up. 
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8 10 



FIG. 3: Variation of 6 (Eq. ^3 as a function of time. The 
parameters are set as follows: 5 = 1 meV, Mhh.hh = Mih,hh = 
0.5 meV, Vxx = 2 meV, e = 0.1, tq = 3.55 ps, ta = 2.55 ps, 
A = 4.5 meV and fi = 8 meV. 

We now move to the second case, where Vxx = 0. 
We use the parameters 5 = 1 meV, Mhh,hh = Mih,hh 
0.5 meV, Vxx = 0, Oq = 8 meV, tq = 4.2 ps, ta = 3 ps, 
A = 3 meV and e = 0.1; in Figs. |3] and we show 
the time variation of 9 and the variation of basis state 
populations. We can see that it is possible to pick up an 
entangling phase (9 = tt) in this case too, and the effect 
now relies solely on the Forster interaction. 



B. Gate fidelity 

The gate fidelity is limited by a number of factors in- 
cluding: corrections to our effective four level QD model, 
spontaneous emission from excited states, coupling of 
the charge carriers to the underlying lattice which is re- 



1.0 
0.9 
0.8 
0.7 

.2 °- 6 
| 0.5 
O0.4 
0.3 
0.2 
0.1 
0.0 





■10 -8 -6 -4 -2 2 4 6 8 1 
t (ps) 







FIG. 4: Variation of the initial state population for each com- 
putational basis state during a gate operation. The parame- 
ters are the same as in Fig. [3 



sponsible for coupling to phonons, and the possibility of 
non-adiabatic transitions due to the finite gate operation 
time. 

Corrections to the effective four level QD model can 
be approximated as exponentially small in the ratio of 
the spectral width of laser pulse to the energy level spac- 
ing. As we discussed in Section III CI this ratio is much 
smaller than unity, and the reduction in the fidelity due 
to excitation to higher energy levels of the QDs can be 
safely ignored. 

It was shown in Ref.|9j that, as well as avoiding difficul- 
ties due to hole mixing, the adiabatic gating scheme also 
avoids further unwanted transitions related to phonon de- 
coherence. However, since the gate is operated in a finite 
time, which one wishes to minimize, there are necessar- 
ily non-adiabatic transitions between the laser dressed 
states. Such non-adiabatic transitions are described by 
the Landau-Zener (LZ) theory^ If we assume a constant 
rate of change for the detuning, i.e. A(t) = At, the 
condition for adiabaticity is given by 2 /A ^> 1, where 
£1 is the energy separation of the levels at closest ap- 
proach. The probability for an unwanted transition by 
P = exp (— 7rfi 2 /4A). There is a simple physical inter- 
pretation of this 9 : r ~ CI/ A is the characteristic time of 
sweep through the resonance, and so the adiabatic con- 
dition is naturally Or S> 1. In the non- linearized version 
of LZ theory the dependence of the unwanted transition 
is still an exponentially small function of both the re- 
ciprocal of the energy difference between the two levels 
(AE a p ~ l/Cl a p where a, (3 are two general QD states) 
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FIG. 5: Variation of 9 (Eq. 14 7^ as a function of time. The 
parameters are set as follows: 5=1 meV, Mhh.hh = Mih,hh = 
0.5 meV, Vxx = 0, f2o = 8 meV, tq — 4.2 ps, ta —S ps, 
Aq =3 meV and e = 0.1. 



FIG. 6: Variation of the initial state population for each com- 
putational basis state during a gate operation. The parame- 
ters are the same as in Fig. [5] 



and the characteristic sweep time r. The proposed adia- 
batic gate scheme is based on the fact that the different 
avoided crossings between the laser dressed states can be 
easily distinguished, i.e. different avoided crossings oc- 
cur for different values of detunings A and energies E 
(see Fig. For the adiabatic gate scheme in the limit 
Mhh.hh = MihMh < Vxx the relevant Rabi frequency 
is the coupling between the states, |00), and |11) when 
A ps 0, which is the biexcitonic shift Vxx- Therefore 
one can estimate the probability of unwanted transitions 
to be around 10 -6 for a typical sweep time of t ~ 4 ps. 
For the case in which the gate is based solely on the 
Forster transfer interaction Mhh.hh = Mih.hh = 0.5 meV 
the fidelity is lower and the probability of unwanted tran- 
sitions for the same typical sweep time is of the order of 
a few percent. In order to get a better gate fidelity for 
the gate based on Forster interaction the typical sweep 
time r should be made longer. 

The adiabatic gate procedure may also reduce the ef- 
fects of phonon decoherenceiS If the gate is operated in 
the bi-excitonic mode the probability of unwanted tran- 
sitions and decrease in pure dephasing decoherence ef- 
fect.s^iiikLL is expected to be of the form 



P 



ft 



cxp(— Aftr) , 



(51) 



where A is a positive constant of order unity and J(u> m ) 
is the spectral function that describes the coupling of 
phonons to our system evaluated at some high-frequency 



cut-off uj m that is imposed by the speed of the frequency 
detuning sweep. 

The typical time scale for spin dephasing is of the order 
of /is and can be made even longer by optical pumping 
of nuclear spins42 Polarizing the nuclear spins will in- 
crease spin coherence since the main mechanism for spin 
dephasing is the coupling of the electronic spin to the 
nuclear spins4^ Thus the ps time scale for the adiabatic 
gate does not constitute too strong a restriction. 



V. STATE PREPARATION, MEASUREMENT 
AND SCALABILITY 

Though our paper focuses on an adiabatic two qubit 
gate scheme which allows one to resolve the difficulties 
arising due to hole mixing, we would like to briefly com- 
ment on the possibility for optically fulfilling the other 
requirements for a quantum information implementation 
scheme, i.e., initial state preparation, measurement, and 
scalability. 

The essential test for any implementation scheme is 
of course provided by experiment. Recently there has 
been tremendous experimental effort and success in val- 
idating the essential stages needed for solid state quan- 
tum computation implementation schemes employing op- 
tically driven charged QDs. In a new experiment Gu- 
rudev Dutt et. al. have demonstrated that a coherent 
optical field can produce coherent electronic spin states 
in QDs, demonstrating that these electronic spin states 
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have life times much longer than the exciton coherence sions. We are indebted to R. G. Beausoleil, who provided 
times. 44 Other experiments have shown how the spin numerical simulation code, 
state of the resident electron in a self-assembled InAs- 
GaAs QD can be written and read using circularly po- 
larized optical pumping. 42,4 ^ Methods for an optical read- 
out mechanism of the spin of an electron confined to a QD 
have also been theoretically suggested^^ 1 ^ and some 
experimental work towards this goal has already been 
performed for colloidal semiconductor QDs.- 49 Scalability 
is a potential problem, since spatial selectively of individ- 
ual qubits is not possible in our system. This is because 
the optical wave length of the exciting laser pulse is much 
larger then the inter-dot distance needed to couple our 
qubits. To optically resolve different QDs we would need 
to resort to energy-selective addressing methods within 
different QD clusters whose size can be controlled^ In- 
side each QD cluster one could spectrally differentiate a 
QD by applying a gate potential and inducing a Stark 
shift of the exciton levels^ Alternatively, globally ap- 
plied pulses could be used within a cellular-automaton 
schemei^^i One other possibility is that small scale pro- 
cessors based on our scheme could be joined together by 
using the exciton coupling to single photons. Linear op- 
tics techniques could then be used to entangle the states 
of two small processors 



VI. SUMMARY 

To summarize, we have derived the Hamiltonian for a 
pair of spin qubits in adjacent QDs, which are coupled by 
trion states. Our analysis shows that it is possible to per- 
form a non-trivial two qubit gate in this system by using 
chirped laser pulses, even in the presence of hole mixing, 
and for two different types of interaction. In Ref. it 
was shown that a dipole-dipole interaction can be used 
to mediate an adiabatic gate of this type. Our extension 
of that approach to include Forster processes general- 
izes this to cover all significant excitonic interactions in 
a coupled dot system. 

In contrast to our previous work^S the gate proposed 
here could be performed in many different materials, 
since we have shown that it is generally valid for all sig- 
nificant excitonic interactions and for varying degrees of 
hole mixing. We therefore believe that demonstration 
experiments could be performed in the near future. 
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